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1. Introduction 

The magnetohydrodynamics of conductive fluids plays an important role in the 
explanation of the existence and stability of the magnetic field of the Earth as well 
as of the magnetic fields of stars and galaxies studied in astrophysics [1] - [15]. In 
the theoretical description of the motion of such fluids an important role is played by 
the induction of magnetic fields which are able to generate a genuine global dynamo 
effect [IB]. For certain spherically symmetric field configurations after a mean field 
approximation a so called a 2 — dynamo model can be obtained [2j. Its description may 
be reduced [2] to the coupled pair of ordinary differential phenomenological equations 

[mill] 

-d 2 r( t>(r)+V u (r)^r)-a(r) X (r)=-X ( p(r), (1) 

- d 2 rX (r) + V d (r) X (r) + d r a{r)d r (j>{r) - V m (r) 0(r) - -A*(r) . (2) 

They are defined on a finite interval of a single coordinate r E (0, R) with R = 1 
after a re-scaling. For the purely kinematic factors one has to set V u (r) — Vd(r) = 
V m (r) / a(r) = £(£ + l)/r 2 where the integer parameter I = 0, 1, . .. coincides with the 
angular mode number of the field. The only input information about the flow of the 
charged fluid (or plasma) is carried by the shape of the so called a— profile a(r). 

Let us now formulate our present task as the construction of the solutions of 
equations |T]) and @ specified by realistic physical boundary conditions 

0(0) = 0, [d r 4>(r)\ r=R + ^(R) = Q, (3) 
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X(Q)= X (R)=0, (4) 

The presentation of our results will start by an outline of consequences of the difference 
between the latter two boundary conditions in section [2l In section [3] we show how 
this difference diminishes in proportion to the quantity \j I. This, of course, indicates 
the possible applicability of a standard perturbation expansion in this parameter. 
After a few introductory technical remarks made on such an approach in the separate 
Appendix A we demonstrate, in the next section [H that the parameter \jl can 
really play the role of a measure of perturbation of the solvable I — > oo limit. In 
the subsequent sections [5] and [5] we then describe the explicit realization of such a 
programme in more detail. Finally, section [7] summarizes all the key ingredients of our 
present new method of solution of the magnetohydrodynamical a 2 — dynamo problem 
by the asymptotic series in terms of certain rational powers of 

2. The dynamo problem in two-space formulation 

There exist several formal analogies of the dynamo problem with quantum mechanics 
[7J GO EH E] • We intend to employ here some of these analogies as a methodical guide. 
A deeper discussion of this aspect of the problem can be found in refs. fTB] - |20j . 
Nevertheless, one should keep in mind that these analogies are violated, first of all, by 
the presence of mixed-type boundary conditions. 

2.1. The doublet of bases 

The most challenging mathematical feature of the physical boundary conditions §3§ 
and ([?]) is that they prescribe the use of different spaces, say, V u and Vd for the 
respective channel functions 0(r) (representing the so called poloidal mode of the 
magnetic field) and x( r ) (representing the complementary, toroidal mode). In a 
compactificd Dirac's bra-ket notation this means that one is forced to employ different 
symbols for elements of each of these spaces. Let us employ the usual vector or ket- 
symbol | x) € Vd in place of the function x(r) emphasizing that for this Hilbert-space 
element the Dirichlet boundary condition Q at r = R = 1 is entirely standard. 
A modified abbreviation | </>} S V u will be introduced for <fi(r) which obeys the less 
standard (mixed, Robin) boundary condition ([3j containing the derivative. 

In the new notation one may introduce orthonormalized bases of eigenvectors 
resulting from the spectral representations of self-adjoint differential operators in any 
of the two vector spaces V u and Vd- We shall postulate 

oo 

- d 2 r + V u (r) = J2 \m}r m {m\, (5) 

oo 

-d 2 r +V d (r) = Y / \n)g n (n\, (6) 

knowing that due to the simplicity of V Ut d(f) and due to the coincidence of the 
boundary conditions in the origin, all the basis functions 4> n (r) = {r\ n} as well as 
Xn( r ) — (r\n) must be proportional to Bessel special functions (ref. [21j may be 
consulted for more details). The basis states have the same closed form with different 
scaling, 

¥>m(r) = M m r 1/2 J e+1/2 (y/r^r) , Xn(r) = N n r 1/2 J f+1/2 (y/g^r) . (7) 
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The kinetic operator eigenvalues r p and g q are given as roots of the appropriate 
combinations of Bessel functions at r = R = 1 . Their numerical values are available 
with arbitrary precision, therefore. 

More care must be paid to the operators which couple the channels. Once they 
represent a map between V u and Vd we can only employ a non-diagonal, full-matrix 
formula for the operator of the multiplication by the function a(r) in eq. ([1]), 

oo 

«(r)= E \k}{k\a(r)\j)(j\. (8) 

k,j=0 

The technique of the evaluation of the matrix elements {k\ a(r)\ j) (by integration) is 
standard. Mutatis mutandis, the same comment applies to the other channel-coupling 
operator 

oo 

d r a{r)d r - V m (r) = - E I k ) ^k,j {j\ (9) 

fej=0 

in eq. ©. 



I!. ,2. Linear algebraic form of the coupled eguations 

The spectral series ([5]) and (|6]) as well as the full-matrix formulae ([8]) and ([9]) should 
be inserted in the coupled set of differential equations ([T]) and ©. Using the pair of 
the natural ansatze 

oo oo 

l<« = El fc }^< lx> = £li>x;> (io) 

k=0 j=0 

we get the set of relations 

oo oo 

E I TO > ( r ™ + A) ro = E I fc i X, , (11) 

m=0 A:,j^0 
oo oo 

E I fc ) Wfc j ^ = E I n ) («» + x )xn. (12) 

fc,J=0 n=0 

Once we assume the completeness of both orthonormalized basis sets { | p] }^L an< i 
{ 1 1) }^Lo spanning respective infinite-dimensional vector spaces V u and Vd we are 
allowed to multiply the previous pair of equations by the respective bra-vectors {p\ 
and (q\ from the left. This leads to the final formulation of the dynamo problem in 
the form of the infinite set of linear algebraic equations, 

oo 

(r p + A) & = £]{p|a(r)|j)xj-, p = %!,..., (13) 

3=0 

OO 

E ^ = (ei + x ) Xq, q = o, 1, ... . (14) 

3=0 

These equations may only be solved numerically (in this case it makes sense to start 
from the elementary elimination, say, of all the unknown quantities <f> p using eq. (|13l) ) 
or perturbatively (in the latter case it seems better to preserve the linearity of the 
whole set in A). 
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2.3. The simplified model with constant a(r) 

Both operators and ^ become perceivably simpler for constant a—profiles. It 
allows us to simplify eq. (fTTj) . 

oc oo 

I k} (r k + A) fe = a I fc } J') K ■ 

k=0 k,j=0 

In parallel, a factorization of the matrix u>kj in eq. (J9j> , 



d r a d r - a £(£+l)/r 2 = 
converts eq. (fT2)) into the simpler relation 



-a \k)(k\j}r 3 {j\ 
k,j=a 



a o E I n ) ( n l m ^ Tm ^ m = E I n ^ ^" + 

n,m— n— 

As a result one gets the pair of coupled linear algebraic equations 

oo 

(r fe + A) 4> k = a ^2 i k \ J) Xj > ( 15 ) 

3=0 

OO 

a 22 (k\ m} T m 4> m = {Qk + A) Xk , fc = 0, 1, ... (16) 

which replaces eqs. (TT3l) and ([M]) at a(r) = ao- For numerical purposes we may easily 
eliminate one of these relations and get the single Feshbach-type set for k = 0, 1, . . ., 



E 



(r fc + A) S k , m - "o E { fc l J') 



1 



j'=o 



or, alternatively, 



E 



«8£{*l. 



m=0 

i.e., after a slight modification, 



3=0 



{Qj + A) 



(r, + A) 



(j|m}< 



(j|m} 



m = (17) 



(18) 



E 

m=0 



(gfc + A - «o) £fe, m + "o E 



A 



fa + A) 



(j|m} 



0. 



We see that the knowledge of the single matrix of overlaps {k\ j) is the only input 
needed for the standard numerical solution of these equations. Unfortunately, the 
latter observation is in fact of no immediate impact upon applications since the choice 
of the constant a(r) — ao ^ is not too realistic (TJ [2]. Even in a more formal 
evaluation it seems oversimplified as it allows many specific relations between its basis 
states 0. One should note that the eigenvalues even become available as roots of a 
certain quadratic superposition of Bcsscl functions as a consequence [2] . In this sense 
our toy problem a(r) = ao remains exactly solvable. Nevertheless, it deserves full 
attention as a very useful methodical guide to more realistic situations. Our present 
ambition will concentrate upon the proposal and description of its new, promising and 
fairly natural perturbative treatment. 
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3. The confluence of bases in spaces V u and Vd at I 3> 1 

There exist several important differences between the present eigenvalue problem (| 13(1 
+ (fl4|) and its quantum- mechanical coupled-channel analogues [22] ■ Firstly, one should 
note the difference in the sign convention (the energies of quantum mechanics would 
be E = — A). Secondly, the present eigenvalue problem is not self-adjoint in the form 
which would be usual in quantum mechanics (consult again ref. 21j for a deeper 
discussion of this aspect). Thirdly, one usually does not encounter Robin boundary 
conditions in quantum mechanics. In this sense our present problem is perceivably 
more complicated. 



3.1. Boundary conditions 

A formal key to the simplification of equations (|13p + (|14p can be sought in a decrease 
of the difference between the Robin and Dirichlct boundary conditions ((3|) and ((4|) at 
r = R = 1 with the growth of £. In a way which significantly weakens the above- 
mentioned constant-profile assumption a(r) = ao let us merely assume now that a(r) 
remains more or less constant in a small vicinity of R = 1. In this case with r « R, 
equations JTJ) and ([2]) acquire the following approximate form, 

-d 2 <f>(r) + (n 2 + A) 4>{r) - a(R) X (r) - 0, k 2 — , 

- d 2 X (r) + (k 2 + A) X (r) + a(R) d 2 r <j)(r) - a{R) k 2 (f>{r) = , r^R. (19) 

This system (with constant coefficients) possesses eigcnfrcqucncics \i = fi[X,K,a(R)) 
obtainable directly from the corresponding characteristic equation 

( M 2 + K 2 + A) 2 = a 2 (i?) {v 2 + k 2 ) . 

In the light of boundary conditions at r = R we may expect that x( r ) = D sin fj,(r — R) 
and <p(r) = C sin (i(r — S) where S ~ R is unknown and where R = 1. An estimate 
of the value of S may be deduced via insertion from eq. ([3]) , 

£ 

HCOs^(R-S) + -smn(R-S) = 0. (20) 
R 

For all I > this enables us to conclude that S > R = 1 so that < Qf. at fc = 0, 1, 

At very large angular excitations the following explicit frequency-independent 
estimate results from eq. (|2U)) . 

5 = l + i +0 (i) 2 . (21) 

We see that the two boundary conditions ([3|) and (j4|) and, hence, also the two bases 
spanning the spaces V u and Vd will coincide in the limit £ — > oo, 

lim r„ = g„, lim | n} — ►|ra), n = 0, 1, . . . . (22) 

« — >oo ^ — ^cx: 

This is an important observation. From the practical point of view it means that we 
shall have 

lim(g|n}=<L „ lim (q\ a(r)\ n} = (q\ a(r)\ n) (23) 

t — *oo I — >oo 

for all the indices q > and n > and for all reasonable functions a(r). 
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3.2. Illustration: Exact solvability of the £ 3> 1 model at a constant a 



As the only dynamical input in eqs. ([TJ) and the a— profile function is a key to all 
the practical and phenomenological considerations. At the same time, all its not too 
large deviations from a constant mean a— profile ao may be treated perturbatively in 
the way outlined in ref. [21] . 

Let us now turn attention to the unperturbed problem with a{r) — a® where we 
add as another assumption that i ^> 1. Due to the limit of the orthogonality rule (f23|) 
the resulting doubly simplified algebraic eqs. p^|) become completely decoupled. At 
every index k the condition of vanishing secular determinant remains trivial, 

A 



det 



-Ot T k 

and its closed solution exists, 



-a 
Qk + A 



0. 



Tk « Qk- 



Ai,2 - ^ 



-Tk 



Qk ± V ( T fe ~~ 8k) 2 + ±Qkdo 



(24) 



In the normalization tfik 

(Xfc)l,2 = T fe + A l,2 



ao the closed form of the coefficients is equally elementary, 
1 



T k - g k ± V [Tk - Qk) + ^Qka? 



(25) 



As long as we have Tk ~ Qk, the latter recipe reproduces exactly the definition 
A = A± = - Qk ± ao^/Qk of the spectrum as obtained earlier under an alternative 
assumption of the high conductivity [H [12] of the dynamo's fluid (cf. eq. (13) in 
ref. [2T]). 



4. Effective simplifications of boundary conditions 

Several complications outlined in section [2] disappear when one returns to the zeroth- 
order I — > oo approximation in eq. (|2ip 21] . For this reason we intend to pay attention 
to the next, first-order level of approximation in l/l. 



4-.1. Approximation using a pair of Dirichlet boundary conditions 

We feel strongly motivated by the observation that eq. (|2ip in fact prescribes an 
"effective" replacement of the Robin (i.e., "difficult") boundary condition ((3|) by the 
Dirichlet (i.e., "easy") boundary condition, say, 

0(0) =0, <f>(S) =0, S = S(£) = l + j. (26) 

One should note that the new condition is specified just by a rightward shift of the end 
of the interval. An important "user-friendly" feature of this shift should be seen in 
the fact that on the first-order level of precision the value of S(£) remains independent 
of all the other parameters. We may expect that the approximation ([3]) — * ([2^]) will 
represent a fairly reliable and nontrivial approximation at all angular mode numbers 
which are not too small. The idea might even remain applicable in the very "realistic" 
domain of the smallest angular mode number £ where a further improvement of the 
choice of S(£) could be sought via eq. (f20|) whenever necessary. 

As next step a (purely formal) transition to all the real axis of r e 1 will be 
mediated by the immersion of both coupled (i.e., poloidal and toroidal) subsystems of 
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the whole system in infinitely deep square wells Vy S ® (r). They are defined as very 
large (or infinite) everywhere except the interval (0, Y) where they should vanish, 

{+oo, r S (— oo, 0), 
0, re(0,K), (27) 

+oo, r £ (Y, oo). 

This means that we change the definition of the kinematic factors in eqs. ([T]) and 
where we set Y = S = 1 + l/£ > I and Y = R= 1, respectively, 

V u (r) = V^\r) = + V^ w \r), 

Vd (r) = V j SQW \r) = + vi 3QW \r) . (28) 

The latter convention is purely formal and its use merely emphasizes the essence of 
the introduction of the approximation (|26p . 

4-2. Models with smeared boundaries 

In the context of paragraph 14.11 the use of the spectral series ^ remains based on 
closed solutions of the Sturm-Liouville problem, 



d 2 r Xn(r) + Xn(r) + V { R SQW \r) Xn {r) = g n X (r) 



(29) 



The trick extends to the Robin boundary condition immediately. The parallel 
construction of the second auxiliary spectral series §5§ is based on the solution of 
the similar equation 

- d 2 Mr) + Mr) + V^ SQW) (r) Mr) = r n 0(r) . (30) 

As long as the effect of the Dirichlet boundary conditions is mimicked by the action 
of an infinitely deep square-well potential, a perceivable reduction of the complexity 
of the basis |n} as well as of the modified space V u is achieved. It is also more easy to 
evaluate quantities t„. 

In summary, the effective simplification of boundary conditions of paragraph 14.11 
looks satisfactory. It is possible to argue that although the numerical performance 
of the resulting simplified bases could be enhanced significantly by the modification 
of boundary conditions, the precision of results may still be kept under control by 
making the modification "infinitesimal" (i.e., controllably small). One can only feel 
dissatisfied by the non-analyticity of square wells ([27]) which does not seem to open 
any easy way towards understanding of the observable simplifications at I — » oo. 

For this reason it would make sense to replace the non-analytic auxiliary square 
wells in eqs. (|29[) and (|30[) by some "infinitesimally" modified analytic approximants. 
One should add that the latter idea need not necessarily be in any conflict even with 
the underlying experimental setup. Moreover, a transition from the discontinuous 
functions Vg S ® W ^ (r) to many of their available respective analytic alternatives 

u[^\r) may be expected better tuned to the analytic essence of perturbation 
methods. 
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In this spirit let us introduce models with smeared boundaries characterized by 
a final redefinition of the corresponding non-analytic functions ()28D . 

V d (r)=vj AA Hr) = £ -^ + U r\r). (31) 
The two requirements of a sufficiently precise fit to the square wells, 

U[ AA) (r) « V*<™ (r) , (r) « vf ^ (r) (32) 

are the only limitations of our free choice of the new analytic auxiliary potentials U. 

5. The idea of l/£ perturbation expansions 

The effective shift of one of the boundary conditions in eq. ()26[) obviously introduces 
a "hidden" kinematical parameter l/£ in the dynamics of our present a 2 — dynamo 
models. The numerical smallness of this parameter is one of our present most 
important observations. We believe that the existence of such a small parameter offers 
an entirely new hint and encouragement for a new development and/or (typically, 
perturbative) simplifications of many existing practical calculations. 

In principle, a prospective way of doing so might consist, e.g., in a systematic 
improvement of the linear-algebraic and weighted-residual approximation techniques 
based on the truncation of the infinite matrices to their finite-size approximants. 
A complementary picture might be provided by the Fourier-series based Galerkin 
techniques of ref. [21] etc. Still, we believe that the best use of the smallness of l/£ in 
our present a 2 -dynamo models can be achieved via a direct use of the various large— £ 
expansions as tested, in the various contexts, e.g., in refs. [23j - [32j . 

Unfortunately, the actual potential of the large— I expansions has only rarely been 
tested out of its natural quantum-mechanical domain. Thus, only this experience 
as outlined briefly in Appendix A below is, at present, available as our preliminary 
methodical guide. Still, we believe that it offers a sufficiently strong encouragement 
of the transfer of the l/£ asymptotic-series techniques to the present MHD eqs. ([T]) 
and ([2|) accompanied by the boundary conditions ([3]) and (J4j> . 

Let us now complement the outline of the method as presented in Appendix A 
by a more detailed account of its features which transcend the routine applications, 
say, of refs. [33). [34] . 

5.1. Coupled channels and the smeared boundary conditions 

In the generic context of X/£— expansions we have to clarify, first of all, the explicit 
perturbation account of the coupling of the poloidal and toroidal modes. In the sequel 
we shall pick up again equations (| 13() and (|14j) with a(r) = ao. Formally, this enables 
us to re-arrange the original differential equations, 

-d^(r)+V u (r)^(r)-a x(r) = -\<f>(r), (33) 

- d 2 rX (r) + V d (r) X (r) - a [a X (r) - A 0(r)] = -A X (r) . (34) 

This simplification will also enhance the transparency of our considerations, paving 
the way towards their extension to more realistic nonconstant a— profiles. 
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In the l/£— series perspective of Appendix A, we have to add now a discussion 
of a few specific features of a 2 dynamos. In a preparatory step let us emphasize that 
the necessary underlying transition to the smeared boundary conditions will imply 
the necessity of the use of the coordinates along the whole real axis, r E K. For 
any smeared-boundary analytic potential (|3T|) we will then be allowed to write the 
corresponding coupled eqs. (|33|) and (|34p in the partitioned operator form 

-d 2 + V u (r) -a 



a X -8* + V d (r) -a 2 ) \ X {r) ) X { X (r) ) ' (35) 

As a consequence, we need not distinguish between the two distinct linear spaces V u 
and V d . We may search for the smeared-boundary solutions corresponding to the low- 
lying eigenvalues —A in the single and standard Hilbert space L2(M). Nevertheless, in 
the perturbation context we shall still be forced to employ the two different bases. 

5.2. Non- coincidence of the two local minima. 

As first step, paralleling the procedure of Appendix A, we make an appropriate choice 
of the two boundary-simulating analytic functions U^ A d A (r) = U Utd {r) of r and expand 

V u , d {r) = ^4r^ + U U:d {r) = V Uid (T U;d ) + V^ d (T u<d ) (r - T Uid ) + .... (36) 



Then the two parallel generalizations of eq. 

K(T„) = 0, 7j(T d ) = 0, 
i.e., the two elementary quadratic equations 

2£(£ + 1) = Tl U' U {T U ) = Tj U' d (T d ) (37) 

define in principle the pair of inverse functions T u>d — T u ^ d (£). By assumption they 
should grow with £, say, in such a way that the two new independent measures of 
smallness 1/T u<d (£) converge to zero in the limit I — > oo. 

The pair of the coordinates T Utd of the respective minima of V Uyd (r) will be 
different in general. In a way paralleling the guidance by Appendix A the respective 
shapes of V U:d (r) near their minima will be specified by the next Taylor-series term 
with the coefficient proportional to the minus-fourth power of the scaling factor, 

T/" (T \ 6 ^ + 1) I TT" (T \ 3 U uA T ^d) , jrll frp \ — rj -4 f^n\ 

V uA T v.,d) = —^4 ^m( J M) = ^ ^ U u, d { Iu > d ) = 2(J u, d - I 38 ) 

1 u, d ±u , d 

In this way we arrive at the zeroth-order descendant of eq. (|35|) , 



a X^a 2 - d i + e +Vd -a 2 cr 2 ) \ x W(r) J { a 2 X ^(r) J (39) 

where we have introduced the additional abbreviation cr^ d V u>d (T Utd ) = v u _ d and where 
we have to keep in mind that the symbol r is just an abbreviation for r = T u + o~ u £, 
or r = T d + a d (. 
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6. Feasibility of perturbation constructions 

6.1. Algebraic form of equations in zeroth-order limit £ — > oo. 

As next step, guided by the notation used in subsection 12.11 we may innovate the 
spectral series ([5]) and ([6]) writing 

oo 

- d l + £ 2 = J2 I m » f ™ m l ' r m = 2m + 1, (40) 



~ d t + ? = Y. \n))Qn((n\, g n = 2n+l. (41) 

n=0 

The specific doubling of the bra-ket symbols has been chosen here to underline the 
specific feature of both sets of eigenvectors: In the zeroth-order perturbation limit 
I — > oo the eigenvalues become well known at all subscripts, Tfc = Qk = 2fc + 1, and 
also the eigenvectors coincide with the two distinct special cases of the harmonic- 
oscillator bases defined in terms of Hermite polynomials [33] . 

It is necessary to add that even in the limit t — » oo the latter two harmonic- 
oscillator basis sets remain distinct in general. A return to Appendix A reveals that 
due to the difference between the Robin and Dirichlet boundary conditions at £ < oo 
the two choices of [say, in the anharmonic form sampled by eq. (fSTj) ] must be 

necessarily non-identical. Hence, also the resulting harmonic-oscillator basis functions 
will differ in both of their ^—dependent shifts T — T u> d(£) and scaling factors o- u ,d- 
This means that we have to re- write our unperturbed eq. 139|) in the spectral-series-likc 
form paralleling eqs. ([8|) and {9j), 



Em=0 I TO }} ( f ™ + w «) ii m \ ~ a l "0 Em,j=0 I 



°d "OA T,n,P=0 I n ))(( n \P}}{{P\ T,n=0 I »)) (fin + ^d ~ ^"o) (("I J V *( r ) 

'(ihr) )' (42) 



Differently said, the resulting matrix problem remains purely numerical even after the 
insertion of the present counterpart 

oo oo 

l<«} = El fc »^> lx»=£li»X. (43) 

fc=0 j=0 

of eqs. (fTDj) . giving the final infinite coupled set of unperturbed equations 

Tk -^ + \) <t>k-aoY,{{k\j))Xj = 0, (44) 
°« / j=0 

+ Xfe + aoA ^ {{k\ m}} <f> m = , fe = 0,l,... (45) 



which is just a smoothed-boundary parallel to eqs. (fT5]) and ([TBI) . 
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6.2. Conditions of solvability of the zeroth-order equations 

For a given pair of smoothed square wells U^f^ (r) which mimic boundary conditions 
the core of the applicability of the large— t perturbation theory may now be identified 
with a guarantee of coincidence of the two bases in the limit I — > oo. Indeed, only 
in such a case one can find the application of the standard textbook algorithms of 
perturbation expansions [M] sufficiently economical and well motivated, especially in 
comparison with the fairly efficient purely numerical methods of linear algebra. 

More explicitly, we require the exact solvability of the zeroth-order version of the 
set of eqs. (|44|) and (|45|) . It may only be achieved when | j}} — > | j)) for all j = 0, 1, . . . 
in the limit t — > oo. In such a case we shall get the orthogonality rule {{fc| j)) — 6k. j so 
that our set of equations (|44| and (|45|) decouples in infinitely many pairs of equations 

-M 

r2 



-~2 A \ "i. - "ii \/, = ( > ■ "ii A i >/, i !- A - - < 

numbered by index k = 0, 1, . . .. With their secular determinant vanishing, 



det( + . 2 ]=0, - I 



wc arrive at the exact solution 



Ai, 2 - ^ 



2 



(46) 



where we have to insert = 2fc + 1 = and a\ = a 2 d . 

Although formulae (|46| look formally very similar to the examples we studied in 
subsection !3.2l the key difference lies in the fact that the coincidence of the bases may 
now be achieved easily. In the light of eq. (|54|) the necessary and sufficient conditions 
of this coincidence | j}} — | j)) reads 

T u = T d = T, V^Tu) = Vf(T d ) . (47) 

Due to eq. ([37)) the first rule means that 

U' U {T) = U' d {T) (48) 

while the validity of eq. (f3"5)) then implies that the rest of eq. (|4"T|) is equivalent to the 
condition 

K(T) = U' d \T) . (49) 

We may conclude that the pair of eqs. (|48)l and (|49|) (i.e., in a geometric language the 
so called osculation of the two curves) represents the necessary and sufficient condition 
of the required coincidence \j}} — | j)) of the zeroth-order bases. 

In practice we may expect that once we fix the value of T 3> 1, we may guarantee 
the solvability of our present illustrative zeroth-order a 2 — dynamo eigenvalue problem 
simply by the choice of the pair of the smooth-boundary-simulating potentials which 
obey the formula 

U u , d (r) = U u<d {T) +A(r-T) + B(r-T) 2 + f u4 {r - T) (r - T) 3 . (50) 

This formula contains the same pair of parameters A and B and two different functions 
fu,v( r — T) which remain regular at r = T as well as two different "offsets" U Utd (T). 
Obviously, the latter formula offers enough freedom for an arbitrarily precise and 
explicit necessary fit (1321) of the boundary conditions. 
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7. Summary 

Despite our explicit knowledge of the exact basis states (0, the immediate linear- 
algebraic recipe (fT3| + (fT4| for the construction of the solutions of differential eqs. ^} 
and © is a numerical task. The necessary matrix elements must be computed as 
integrals of the products of the pairs of Bessel functions with given and, in principle, 
arbitrary phenomcnological input function a(r). The matrix elements rarely remain 
available in closed form and sophisticated numerical methods are necessary for their 
evaluation. Moreover, even if the precision of the matrix elements themselves proves 
satisfactory for a given a(r) in ((5|) and ©, the final solution of the eigenvalue problem 
requires an infinite-dimensional matrix diagonalization. 

In this context we revealed that the quantity l/£ represents a "hidden" natural 
small parameter in the problem. This encouraged us to search for an efficient non- 
numerical construction of the solutions. An ambitious candidate has been sought and 
found in the perturbation method of the so called large— t expansions. A detailed 
adaptation of the key ingredients of this technique to the specific needs of the coupled 
differential a 2 — dynamo equations has been performed here in some detail. 

As a first step of the realization of such a project the purely Dirichlet specification 
of the toroidal mode x( r ) has been selected and interpreted as if resulting from the 
action of an infinitely deep square well V^ s< ^ w '(r) which vanishes precisely inside the 
corresponding finite interval of r. This trick gives an equivalent picture and at the same 
time it offers a guide to the simplification of the treatment of the more complicated 
poloidal mode 4>(r) constrained by the Robin boundary condition. In this sense the 
Robin constraint has been re-interpreted as the same Dirichlet boundary condition 
modified by its shift to some slightly higher value of an "effective" boundary point. 

In the second step towards the realization of perturbation solutions a "softening" 
of the "rigid" boundary conditions has been introduced via a replacement of the "deep 
well" V( S Q W ' (r) by an element of a large family of its suitable analytic analogues and 
descendants V^ AA '{r). For the sake of definitcncss the effects and consequences of this 
"softening" have been illustrated on the power- law well V {AA \r) ~ £ ax x with a 
suitable, not necessarily integer exponent K 3> 1. We reminded the readers how the 
1/i technique works in this older and, by the way, numerically extremely successful 
[26] application. 

In the last third of our paper we demonstrated that and how the "softening" of 
the boundary conditions opens a mathematically consistent path towards the use of 
the inverse mode number I as the effective small parameter of the theory of spherically 
symmetric a 2 — dynamos. In particular we showed that the asymmetric character of 
the coupling of the toroidal and poloidal channels leads to certain unexpected and 
nontrivial formal challenges and we were successful in finding their resolution in the 
use of a certain flexibility in the choice of the softened square wells. 

Let us add in conclusion that the details and technical aspects of our present 
method of construction will strongly depend on assumptions concerning the structure 
of the input a— profile a(r). In fact, a strong sensitivity of the results on such 
details has been observed and studied using a more standard perturbation theory in 
ref. [21]. Our present text skipped the details of all the wealth of phenomena related 
to the variations of a(r). In the light of our present emphasis on the description and 
clarification of new methods we often restricted our attention to the mere constant 
choice of a(r) — ao. This limitation seems to have left a lot of space for its 
phcnomenologically motivated weakening or even complete removal. 
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Appendix A: l/£ expansions in Quantum Mechanics 

Quantum anharmonic oscillators and their Schrodinger equations 

[-8 2 r + V(r)] ip{r) =Ql>(r), V(r) = W+H + ^ r 2K (51) 

are often studied in the domain of I ^> 1. Reviews [26] and [27] of the corresponding 
techniques may be consulted for many technical details. In a sketchy outline of these 
techniques let us emphasize that eq. (|51|) might be also reinterpreted, in our present 
a 2 — dynamo context, as one of the eligible smooth-boundary alternatives to eq. ([25]) 
or ([30]) . The single-channel Sturm-Liouville problem ([ST]) defines then the basis states 
\m) or \n} entering the spectral series ([6]) or ([5]), respectively. 

In a preparatory step of the current 1 /£ recipe let us recollect that our function 
V(r) may be re-written as the Taylor series 

= g(£ + l) + u 2 r K = y(T) + V'(T) (r-T) + l V"(T) (r — T) 2 + . . . . (52) 
r z 2 

At a point T defined as the unique global minimum of V(r) we have, in particular, 

V'(T)=0 =j. 2 ^t 1} ^Koj 2 ^- 1 =j. T k+2 = 2 ^ +1) ■ (53) 

This means that for large i the coordinate T = T(£) is also large. We may use the new 
measure of smallness 1 /T in place of the original small parameter 1 / £ in the formulae 

Ku 2 T K+2 

V(T) = — + ui 2 T K = lo 2 T K (1 + K/2) , V'(T) = , 

11 z 

V"(T) = Klj 2 T k - 2 (K + 2) , V"'(T) = Klj 2 T k - 3 (K 2 -3K-10), .... 
We insert them in eqs. ([52]) and (]51[) and change variables in such a way that 

r " T = f ^' ff4 = k(k + 2 5 Ur*-2 • ( 54 ) 

This reduces eq. (|5ip to the formally equivalent equation 

-fl^(r) + K{K 2 + 2) oj 2 T k ~ 2 (r - T) 2 ^(r)+ 

+ X(K2 ~!^~ 10) - 2 ^- 3 (r - T) 3 ^(r) + ... = [,- V(T)] ^(r) (55) 
6 

and, after rescaling, to the perturbed harmonic oscillator, 

[-d 2 + e + K 3 e + K4ti 4 + ...] i/>(T + irZ)=e1>(T + tTS). (56) 

In the zeroth-order approximation its low-lying unperturbed spectrum is well known 
and equidistant, 

e = a 2 [£>-F(T)]«£ = l,3,5,... 

and all the higher-order perturbation corrections may be evaluated easily |33j . 
Moreover, their size may be made arbitrarily small by the choice of a sufficiently 
large I since it is easy to show that we have n 3 — O (T~ 1 / 2 ~ K / 4 ) , k 4 = O {T~ 1 ~ K / 2 ) 
etc. This means that we get the pure and exactly solvable harmonic oscillator in the 
limit £ — > oo. Perturbation theory helps us then to evaluate the corrections at any 
finite £ < oo and for all the low-lying levels with n = 0, 1, . . ., 

Qn = V{T) + a- 2 (e„ +£! + ...) (57) 
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with Ej = 0(K j+2 ) = O (T-3/ 2 -J K / 4 ) for all perturbation orders j = 1,2,..., i.e., 
with 

g n = lu 2 (1 + K/2) T K + uj [2n + 1) \J K(l + K/2) t k / 2 ^ 1 + ... 

in our illustrative example. 

Marginally, let us note that the fact of the large size of the shift of the position 
of the global minimum T > 1 of the complete, so called "effective" potential term 
V e ff(r) = i{l + l)/r 2 + V{r) in the Schrodinger equation resolves, as a byproduct, 
also the well known puzzle that our approximative wave functions are, in general, 
allowed to become singular near the pole of the centrifugal component of V e ff(r), i.e., 
near r = 0. In the present approximation-theory context, the practical irrelevance of 
the corresponding error terms is an interesting consequence of the fact that all our 
approximate wave functions are exponentially decaying outside the effective-potential 
valley at the sufficiently large I. This makes the irregularity of the wave functions at 
r = -C T (which is definitely there of course, together with the equally puzzling 
"allowed tunnelling" to r < [3D]) entirely irrelevant. Still, it is worth noting 
that all these "omissions" definitely cause the ultimate divergence of the infinite l/l 
expansions, in the manner discussed thoroughly in the specialized literature (of which 
we may recommend the most succinct study [26j to the interested reader's attention). 



